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04 ■ Abstract. Let C be the transfer operator associated with the Gauss' continued fraction map, known 

^ — ^ I also as the Gauss-Kuzmm-Wirsing operator, acting on the Banach space. In this work we prove an 

■ asymptotic formula for the eigenvalues of £. This settles, in a stronger form, the conjectures of D. Mayer 
_ " and G. Roepstorff (1988), A.J. MacLeod (1992), Ph. Flajolet and B. Vallee (1995), also supported 

' by several other authors. Further, we find an exact series for the eigenvalues, which also gives the 

■ canonical decomposition of trace formulas due to D. Mayer (1976) and K.L Babenko (1978). This 
00 , crystallizes the contribution of each individual eigenvalue in the trace formulas. 
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1. Conjectures and the main result 



I Our first attack of this problem, based on tlie techniques developed in [1], is presented in [2]. Though 

J3 ' the obtained result is highly supported by numerical computations, the series for the eigenvalues 
in [2] is too complex and does not give neither the structure, nor asymptotics. Here we present 
another approach which does answer these questions. In the second part of this paper we explicitely 
^ ■ compute the first constants in the asymptotic expansion and investigate a refined arithmetic structure 



of functions W£(X). 



. 1.1. Introduction. Let D be the disc {x G C : |z— 1| < |}. Let V be the Banach space of 
. functions which are analytic in D and are continuous in its closure, with the supremum norm. The 
Perron- Frobenius, or the transfer operator for the Gauss' continued fraction map, also called the 
Gauss-Kuzmin-Wirsing operator, is defined for functions / G V by [11, 12, 13, 18, 24] 



^ 1 „/ 1 
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. uautis-ix uzirmi- vvirs'ing open 

r.> ■ m=l 
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5r ■ Our chief interest is the point spectrum of this operator. As was shown in [3, 4], this operator is of 
trace class and is nuclear of order 0. Thus, it possesses the eigenvalues A^, n G N, which are real 
numbers, |An| > |A„+i|, Ai = 1, and Yl'^^=i l-^nl*^ < every e > 0. In fact, the Hilbert-Schmidt 

operator /C, defined by 

00 

^umx)= [ —Jl2^^=u{y)dy, 
^ ^'^"^ ' J ^(e-_i)(ey-l) 

for u belonging to the Hilbert space L'^{M.^,m), dm{y) = dy, has the same point spectrum, 
counting algebraic multiplicities [3, 17, 18]. The tranfer operator C arises from the Gauss map F{x) = 
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{1/x}, X G (0, 1], F(0) = 0. Let F(^) = F, and F^^) = FoF'^^^^^ for A: > 2. As is now well-known, due 
to important contributions by Gauss, Kuzmin, Levy, Wirsing, Babenko, Babenko and Jur'ev, Mayer, 
we have 

M(a e [0, 1] : F^^\a) < z) = ^^-^ + £ \i^n{z). 

^ n=2 

Here stands for the Lebesgue measure, and for each n > 2, the eigenfunction ^n{z) is defined 
in the cut plane C \ (— oo, —1], it satisfies the boundary conditions ^>n(0) = ^n(l) = 0, the regularity 
condition 

sup \{z + l)U{z)\ < +00. (2) 

where U{z) = ^n{z), and the functional equation 

+ = i-.c|>„,(^). 

Thus, ^i{z) = —^^Y^- The eigenfunctions of C are then given by ^'^(z), n G N. Moreover, every 
function U{z) which satisfies, for a certain A G M \ {0}, the functional equation 

U{z) = U{z + l) + j^^^u[j^), zgC\ (-00,-1], (3) 

and the regularity property (2), is the egenfunction of C with the eigenvalue A. More details can be 
found in [3, 13, 24]. 



The nature of the eigenvalues A„ for n > 2 is unknown. It is widely believed that these constants 
are unrelated to other most important constants in mathematics; in particular, it is expected that 
they are neither algebraic numbers nor periods. Now, more that 480 digits of A2 have been calculated 
[5], but one can get rigorous certificates only for the several first digits of A2 and A3 [8, 13, 15, 17, 25]. 
On the other hand, the trace of the operator C can be given explicitly. As was shown in [16] (see also 
[6, 7, 14, 17, 18]), we have 

00 

-(^) - |:^-/|^^-E^4-^-i|(-)'-'(T)<c(-)-). 



here 



^i = 7-JTT, ' = ' ^;^'ieN. 

£+£+£+ i+J+i+J+ 

Without going into detail, we note that these trace formulas and this field is deeply and intricately 
ralated to the Selberg zeta function, the Riemann zeta function, Maass wave forms and modular forms 
for the full modular group [14, 19]. 
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Throughout this paper, we fix the notation 

'^=^— 

The constants have received a considerable amount of attention in recent decades. Nevertheless, 
there were three outstanding unresolved problems; no theoretical progress was made towards any of 
them. As was said before, we henceforth arrange the eigenvalues according to their absolute value 
I All > IA2I > • • • . Of course, in case A„ = ibA„+i for some n, this arrangement is not uniquely defined. 
Despite this, we have 

Conjecture. The following three statements are true: 

i) Simplicity. The eigenvalues are simple, |A„| strictly decreases. 

ii) Sign. The eigenvalues have alternating sign: (— 1)"~^"'^A„, > 0. 
in) Ratio. There exists lim = -^^^ = -4>^- 

The first and the second conjectures can be attributed to Mayer and Roepstorff [18], also reiterated 
by MacLeod, and the last was raised by MacLeod [15] (only with a constant ~ —2.6), and seconded 
by Flajolet and Vallee [6, 8, 9, 10] (with the constant —0^). Several authors also claimed to believe 
these conjectures. The ratio conjecture has the following explanation. The spectrum of the operator 

£o[/WIM = ^/(;^), /.V. 

is given by (— 1)"'+^(/)~^", n E N. It is expected that the terms in (1) for m > 2 act only as small 
perturbations to Cq. Of course, the "Sign" and "Simplicity" conjectures follow from the "Ratio" 
conjecture for sufficiently large n (provided it is effective and we can verify these conjectures for the 
first values of n). 

1.2. Main results. It is surprising that in fact the real asymptotics of the sequence A„, minding the 
above remark about Cq, is as simple as is allowed. 

Theorem 1 (Asymptotics). We have the formula 

(-l)"+U„ = (^-2" + c(n) 



where 0.4 < c{n) < 1.7, and c{n) tends to a limit, as n — )• 00. Further, there exists an asymptotic 
expansion of the form 



00 



{-iT+^Xn ~ <^-'" + d{p)n-P^\ d{p) E M. 

Based on the computation of P. Sebah [22] , we have: 

c(l) = 1.618+, c(2) = 1.529+, c(3) = 1.403+, c(10) = 1.223+, c(20) = 1.184+, 
c(50) = 1.153+, c(70) = 1.145+, c(lOO) = 1.137+, c(149) = 1.1313+, c(150) = 1.1312+. 

Our second result gives the structure of the eigenvalues. Let Pm'^\x) stand for the classical Jacobi 
polynomials [23]. 
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Theorem 2 (Arithmetic and decomposition of trace formulas). There exist functions W„(X), v > 0, 
defined by M^(X) = 1, M^i(X) = | • (^~^-'^i-*^^|(3/2), and then by a certain explicit recurrence - this 
will be given later, see (7) - such that 

oo „ 
1=0 

for an absolute constant C. This decomposition is compatible and gives the decomposition of trace 
formulas for the powers of C: for the first and second powers, we have, respectively, 

oo 

n=l '^i + ^ 

oo -. 

Analogously for higher powers of C. 



So, in the trace formulas now we are able to crystallize the contribution of each individual eigenvalue. 
Thus, this defines an infinite matrix whose elements in rows add up to eigenvalues, and elements in 
columns add up to (^^^ + 1)"^. The last two equalities can be checked to hold true for, say, £ = 1,2 
(in the first case) and £ = 3 (in the second), if we use the representation of Jacobi polynomials in 
terms of integral (6): 

oo 



n=l ^J- 



oo 



MAPLE re-confirms this, too. Further, the last theorem, for the first time, gives rigorous certificates 
(though considerable computational time and space resources are needed) to calculate numerically the 
first few digits of an eigenvalue with any index. 

Corollary. All three claims of the Conjecture are true. 

2. Preliminaries 

2.1. Jacobi polynomials. For a,/3 G Z, m G No, the classical Jacobi polynomials are given by [23] 

here a small contour winds w = x in the positive direction. For arbitrary x outside the closed interval 
[—1, 1], one has an asymptotic formula 



p(«,/3)r , ((x+i)^/^ + (x-i)v^)"+" 1 / 2_..i/2r+^/^ 

™ {X- l)"/2(x + l)/3/2(x2 - 1)1/4 • ^/2^ -[X + ^X I) J 



as m — 7- oo. 
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We will now need the case m = X — 1, (a,/3) = (0, 1) (see Theorem 2). Moreover, one can extract 
the second asymptotic term. In this particular case this reads as 

-Px-i(3/2)~^-^ + 0^^ 

(Here and in the sequal one can think of X = n, exactly the index of an eigenvalue A^, though we 
choose to use an unspecified variable to denote that it is a function in X). 

2.2. g-coefRcients of the analytic function. We will now introduce some special coefficients. 

Proposition 1. Every analytic in the half-plane 5i(z) > — ^ function f{z) can be expanded in the 
following way: 

f^'^ = 2^^-^ (, + ^).+i ' 

where \aj\ < C{f, e) • (1 + e)-^ for every e > 0. We call aj the jth golden coefficient, or (7-coefficient, of 
the analytic function f{z). 

Proof. Note a simple identity 

c ^ 
here a small contour C winds the point cj)^"^ in the positive direcion. So, we define jth (^-coefficient of 
the function f{z) by the formula 



^ 2^ J ^ ' {z- (t)^^y 

Let 



z — cp wcp + 

w 



Z + (j) —w + 1 

The z-half plane 3?(z) > — ^ is mapped in the -w-plane to the disc \w\ < 1. Given /, analytic in 
3^(z) > -i. Let 

\ w — 1 / \ —w + 1 / 

Then g{w) is analytic inside < 1 and Oj is its Taylor coefficient at w^~^. □ 

For the example, the ^-coefficients of the dominant eigenfunction of C, namely, f{z) = are 
given by 

2(t)-l f 1 {z + </>)J 



27ri J z + l (z-0^iy 



dz. 



Expand the contour to a large circle. The residue of the function under the integral at z = — 1 is equal 
to (-1)^V>-2J. Make a substitution z = ^ - I. We are left to calculate the residue aX w = Q. This 
gives, for j > 1, 

Oj = {2(t) - 1) • (1 - (-l)^V)-2i) ~ 2(/) - 1, as j ^ oo. 
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We will need the following crucial result. 

Proposition 2. Let ^ G N. The jth g-coefficient of the function 

{z + l + <pf+^ 

is given by 

K(, ^ 5 . 2^-'-' . r'-' . (3/2) ^ 1^ / d.; (.) 

the contour winds 3/2 in the positive direction. Moreover, we have the symmetry 

iK{j,£)=jK{e,j), j,e>i. 

These coefficients are positive, and ^^=iK{jA) = 1 for every G N. Further, for fixed i, K{j,i) 
achieves its maximum at j = I, and 

5I/4 

K(i,i) -. 

Further, K{j,i) is of fast decay when \ j — £\ increases} 
Proof. By the above remark, 

c 

here C is a small circle around in the positive direction. Let a, b, c, d be distinct real numbers. We 
will generally explore the integral 

nz + ayjz + by-'^ dz 
J {z + cy{z + dy-^ {z + df' 

here the contour winds — c in the positive direction. Let 

z + h z + a w{d — a) + {a — b) p w + q 
= w ^ = — — r 7 — = - • . 

z + d z + c w[a — c) + [c — b) r w + s 

So, after this change, the integral transforms into 

p' f w^~^(w + qy , 
^ = i —r^ dw. 

ri{d-b) J {w + sy 

Further, let us make the change w 1— t- {q/2)w — q/2. This gives 



r^d-b) J {w- l+2s/qy 

In our case, a = (j), b = 1 — 0"^, c = —0"^, d = 1 + 4>. So, 

a — 6_iC — 61 3 

p = d-a = l, r = d - c = 2(j), q = - = 20 ,s = - = -7r[, -l + 2s/g = --. 

d — a d — c 2(j) 2 

So, 

'^'^ I {w-iy-\w + iy ^ 

dw. 



2^(20-1)7 (u;-3/2)J 



^The needed technical details about K{j,i) will be filled in the next version of this manuscript. 
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where the contour goes around the point w = 3/2 in the positive direction. This integral can be 
expressed in terms of Jacobi polynomials, as the definition (4) shows. In our case, m = j — 1, 
a = i — j,l3 = l, x = 3/2, and this gives the statement of the proposition. □ 

3. The proof 

3.1. The main recurrence. As before, let X G N. Let us define 

v<«)(x) = i, „(»)(x,.) = (£^l^, 

and then the functions ^('"^(X), u^'^\'K,z) recurrently by 

^ f 1 \ X+l J.2 X 1 

Y: n'^-'^ ( X, z) v(^) ( X) = ^^f-r^u^"^ ( X' ttt) + ^ ""^^'"'"'^ ( X, . + i)y (^) ( X) . (7) 

■r=0 "T J ^ + j,=o 

This recursion is the core idea of this work: for the unknown function u^^^ ( X, z) we have two terms 
rather than three, since the last sum is delayed. Also, it turns out that V^^^ is the dominant contrib- 
utor in the asymptotics of A„. This recursion works as follows. If we have already determined V^"^^ 
and u^'^^ for r < v — 1, the above identity allows to uniquely determine (7-coefficients of u^'"\ except 
the Xth coefficient, which can be defined arbitrarily; we set it to be equal to 0. This recursion also 
yields the unique value for V^'"'^ , as we will soon see. 

To see clearer the main idea of this paper, multiply the identity (7) by and sum over 

u > 0. If we put 

oo oo 

A( X) = J] </>-2 . y W ( X) , [/( X, = </.-2 x^nW ( X, z) , 

£=0 1=0 

we obtain the avatar of the initial functional equation 

(_l)X+i I 
A(X);7(X,z) = A(X)[/(X,z + l)+ ^. [/(X, 



{z + lf \ ' z + 1 
We will soon show that 

Proposition 3. If aj, j > 1, are the g-coefficients of f{z), then the g- coefficients of 

(_1)X+1^2X 1 



(Z + 1)2 + 

are given by bj = aj{—l)-^~^^4>'^-^~'^^. The g-coefficients of f{z + 1) are given by Cj = Yli^i "i-^(j>^)- 

Let the g-coefficients of u*^^)(X,z) be given by J > 1 (we ommit indication of the dependancy 
on X). Now, let us compare the jth (^-coefficient of (7). We obtain: 

V v—1 oo 

Y: ^r^V('') = qf (-1) x+.^2X-2, + ^ q^-^-'^Kij, i)V^^\ (8) 

r=0 r=0 1=1 

For example, when f = 0, this reads as 

„(0)y(0)=g{.0)(_l)X+.^2X-2,^ 
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and we readily obtain that only is non-zero; in our case, it is equal to 1. When v = 1, (8) reads as 

+ ,(0)y(l) = ,(.l)(-l)X-l-.^2X-2, +^2X^^(^. X). (9) 

When j = X, we obtain 

= (/,2Xk(X, X). 

(Note that this will be the second largest contributor to the value of A^, and it is equal to cl)~'^^K{n, n), 
which is of size - this guarantees the success of our approach!). When j / X, (9) gives 

(1) _ cP^'^Kij, X) 



1 - (-l)X+i,/,2X-2i' 



As we now see, can be defined arbitrarily. Indeed, choosing another value for £ > 1, leads 



to a function C/( X, z) which is diferrent from C/( X, z) by a constant factor X^£>o 9x P^^ovided 
the last series is absolutely convergent. We therefore always choose q-^ = for i > 1. 
As a next step, let us read the recurrence (8) in case v = 2. This gives: 

oo 

1=1 

When j = X, this yields 

y(^) = /-i.^(X, X) + 0- j: r^^J'^'^ 



('_iy+X^2X-2i ■ 
j^X ^ ^ 



(2) 

Choosing j 7^ X gives the value for q^ : 



(2) ^ ,4X X)K(X, X) 

(1 _ (_l)i+X02j-2X)(l _ (_l)i+X<^2X-2j) 

+ ^4XX- K{j,i)K{i,X 



5. fl - (-1) 



(1 - (-l)X+i02X-2i)n _ (_l)X+i02X-2i)- 

X 

The recur rsion works the same way for v > 3. 

3.2. Asymptotics and convergence. ^ We write /(X) =5((X), if 

|/(X)-5(X)| <^ 

for some universal constant C. Let us take a look at (8) and consider this as equality in the sense of 
=. Thus, we care only about terms which are of size By induction, we obtain 



00 



qf)+qf)v(^) ^ g('^)(-l)X+.^2X-2, ^^2X^^{.-1)^(^. ^) ^^2X^(^. x)f(-i). (10) 



i=l 

As we see from Proposition 2, for v > 2 the last term drops down, too. When j = X, this then gives 

00 

(1,-1) 

i=l 



y(-)-02X^^(.-i)^(X,i), v>2. 



^The needed technical details for this subsection - about the asymptotics of K{i, £) for example - will be filled in the 
next version of this manuscript. 
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Next, (10) for J 7^ X gives 

l_(_l)XT3<^2X-2i 



These equalites and the asymptotic properties of the coefficients K{i, 1) (Proposition 2) imply that 
|y(^)(X)| < . ^ ^2X. . C 

for a certain universal C. 

So, we are only left with the following question: how we can be sure that (-l)"+^A(n) gives all 
eigenvalues of £, that there are no sporadic ones? The answer is provided by the trace formulas, and 
for this purpose we introduce a generalized operator. 

3.3. Generalized operator. Let w € C, | < 1. Now, let us consider the following linear operator 
/ : V 1-^ V, defined by 

oo m— 1 1 

^M{mz) = Y.r^—^f{-— 

m=l ^ ' 

This operator is a generalization of (1): Ci[f{t)]{z) = C[f(t)]{z). If G{u),z) is the eigenfunction of 
this operator with the eigenvalue A(t<j), then it satisfies the regularity condition (2), and the functional 
equation 

X{u)G{u},z) = uXiuj)Giuj,z + l) + -^-—-G(u;,^—), z G C \ (-oo, -1]. 

(z + 1)^ V z + 1 / 

Multiply now the recurrence (7) by lo^. We see that, if we put 



£=0 

then A„(a;) = (— 1)"+-'^A( n) is the eigenvalue of C^} with the eigenfunction U{n, u},z), where 

oo 

UiX, a;,z) = J]</.-2X^.^^.^W(X,^). 



Similarly like £i, for fixed cj, \uj\ < 1, the operator is of trace class and is nuclear of order zero 
[16]. Thus, 



°o , ,m.-l °° , ,i+j-2 



5J = 

Now, we know that the eigenvalues of are analytic functions of oj for \u}\ < 1, and they are real 
for real uj. Let the set of eigenvalues be the union of {Xn{uj),n £ N}, and {ai{uj),i £ J^}, where 
is a finite or a countable set. We have: 

oo 

J]A2(a;) + J]a2(a;) = Tr(£2j, |^|<i. 

n=l tej^ 
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In particular, for a; = 0, this gives 

I oo 

E-^'(o) = ^-E^n(o) = o. 

ie,y n=l 

Thus, since o"j(0) G M are eigenvalues of Cq, this implies = 0. Also, comparing the coefiiecients at 
the powers of lo in the trace formula, we obtain 

oo ^ 

5](-irV"^^"^(^"^)(n) = -^, i>l. 

n=l V + ^ 

gives the penultimate identity of the Theorem 2. Also, if we compare 
the coefficient at uj^'"^ at the trace of the operator we obtain 

oo ^ 
i+j=e n=l i+j=£ ^^^'3^:>'^' 
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